Two-photon exchange at low Q 2 



Dmitry Borisyuk and Alexander KobushkirQ 
Bogolyubov Institute for Theoretical Physics 
Metrologicheskaya street H-B, 03143, Kiev, Ukraine 
(Dated: February 4, 2008) 

We study two-photon exchange for elastic electron-proton scattering at low Q 2 . Compact approx- 
imate formulae for the amplitudes are obtained. Numerical calculations are done for Q 2 < 0.1 GeV 2 
with several realistic form factor parameterizations, yielding similar results. They indicate that the 
corrections to magnetic form factor can visibly affect cross-section and proton radii. For low-Q 2 
electron-neutron scattering two-photon exchange corrections are shown to be negligibly small. 



During last few years, two-photon exchange (TPE) in clastic clcctron-proton scattering draw a lot of attention. 
The study of TPE was inspired by the problem in the proton form factor (FF) measurements However, this is 
not the only field where TPE effects can manifest oneself. The others include single-spin asymmetries @, radiative 
corrections to parity-violating observables [|[ and determination of proton radii [j, Q . The latter is done via study 
of elastic electron-proton scattering at low Q 2 . The TPE corrections are known to be important here, but the whole 
TPE amplitude was never calculated in closed analytic form. 

The TPE amplitude can be split into elastic (with proton intermediate state) and inelastic contributions. Naturally, 
the elastic part should be dominant at low energy. In Ref . [|[ the elastic contribution was reduced to the twofold integral 
*vO ' containing FFs from the space-like region only. Nevertheless, the functions under the integral are complicated and 
possess singularities, thus it is not so easy to integrate numerically. At low energies, further simplification is possible, 
yielding compact analytic expressions for the TPE amplitudes. 

We follow the notations of Ref. Q. The particle momenta are defined according to e(k) + p(p) — > e(fc') + p(p'). Wc 
also introduce P = (p + p')/2, K = (k + k')/2 and q = p' — p. Scalar kinematic variables are v = APK = s — u and 
t = q 2 = -Q 2 . 

Neglecting the electron mass, the scattering amplitude can be written as 
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The invariant amplitudes Fi are functions of v and t. In the Born approximation (BA) Fi(v,t) — F\(t), F2{v,t) = 
\ F2(t), F^,(u,t) = 0, where F\,2{t) arc usual Dirac and Pauli proton FFs. Alternatively, wc may use any three 
independent linear combinations of J^; e.g. instead of F\ and F 2 one may define Gm = Fi+F 2 and Ge = F\ + j^F 2 , 
^ . such that in BA Ge and Gm are equal to electric and magnetic FFs. 

The difference between F^ and their Born values is proportional to a ~ y^y. Neglecting the terms of order a 2 w.r.t. 
the leading one, the cross-section for unpolarized particles is 
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da = da (sG 2 E - jhtGlt) , (2) 



Gm + £ 4M? F3 > Ge — Ge + jjp-F 3 , (3) 



8rra 2 AM 2 dt , , 

d(T0= l-eJ^WT (4) 

The invariant amplitudes are complex; in Eq.@ and everywhere below their real parts are understood. 

Eq.([2]) looks like Rosenbluth formula (certainly the experimental separation of Gm and Ge is not possible, since 
both are e-dependent beyond BA). Because of this simple form, we adopt the following set of invariant amplitudes: 
Ge, Gm, and G M . 

In this paper we are interesting in TPE contributions SG = {SGe,^Gm,^Gm) at low t, —t <C 4M 2 . On the other 
hand, we neglect the electron mass m, so it should be Am 2 <C —t. 
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Low t can be achieved in two ways: (i) at fixed energy and small scattering angle and (ii) at fixed angle (or fixed 
e, which is nearly the same) and low energy. We will consider the latter case. This implies 
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v = xJ-t(AM 2 - t) ~ where x=J ——. (5) 

V 1 — e 

In the Breit frame x = 1/ sin | ; this quantity was denoted x in Ref . @] ■ Eq. © suggests that at low t one should 
primarily consider the amplitude Qe- However for generality we will study all three amplitudes. The role of amplitude 
Qm increases for backward angles (e w 0) and for neutron target (Ge ~ 0). The amplitude Gm may contribute to 
polarization observables. 

The TPE contribution to any of the invariant amplitudes can be written as @ 



ss = 5* E/ E ^(tiMm^ih)^^, (6) 

where n = 1, 2, 3, 4, 4x; 



Di = l, ^ 2 = (P-p" + A') 2 -m 2 , £> 3 =p" 2 -M 2 , (7) 
£ 4 = [(P -p" + A) 2 - m 2 ][p" 2 - M 2 ], £ 4x = [(P -p" - Kf - m 2 ][p" 2 - M 2 ] 

and ti = (p — p") 2 , t 2 = ip' —p") 2 . The factor j% is introduced for convenience. A n ^j are some polynomials in 
ti,t 2 , depending also on v and t, and originate from the numerators of TPE diagrams. They are different for different 
amplitudes in the l.h.s., but to simplify notation we do not attach an index to indicate this. The crossing symmetry 
implies 

4i,2,3(i/) = -A 1Aa (-v), Aa{v) = -A ix (-v). (8) 
Following Ref.@, the integrals over d 4 p" can be reduced to the twofold integrals over t\, t 2 : 
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where K, n are known functions, given in Eqs. (31-35) of Ref.[||. It is convenient to define A±± = A4 ± A 4x , Ka± = 
^(/C 4 ± /C 4x ), then 

AilCi + A 4x /C 4x = A i+ K. 4+ + A 4 _/C 4 _. (10) 

Now we will simplify all the terms of Eq.@ in the case of low t. Due to 8- functions, contained in ZC„, the integration 
in ([9]) is done over three regions: (I) Xoo > 0, (II) x^ < < %, t\ + t 2 — t > and (III) xm < < x m , ti + t 2 — t > 0, 
see Fig. [TJ In the regions II and III tip ~ t but in the region I ti and t 2 can be large. However the FFs entering (J9j) 
decrease rapidly as tip — * —00, so the main contribution in the region I also comes from tip ~ t. Assuming tip ~ i, 
we obtain 

JCi ~ 1, £ 2 ~ /C 4 _ - l/i\/ = i, (11) 



AC 3 w — — =[0(0^) + 20(-x oo )0(x M )0(ii+t 2 -i)] « — r^ Oixoc) ~ 1/V=t, 
4m y — t 4M v — t 

2 ■ 2 ■ 

IC i+ « ^^sign(t 1 +t 2 -i)[0(a; oo ) + 20(-a; oo )0(a;M)e(ti+t2-i)] « -^0(a-oo)sign(i 1 +i 2 -t) - l/tv 7 ^, 

where v 2 ( tl+ * 2 ~* ) 2 — <it 2 (^ 2 + 4M 2 t), and the square root means arithmetic value. Since the area of region II 
is small as 0(t), the term 9{— x oo )9(xM)0iti+t 2 — t) brings, after the integration, an additional factor oft. Therefore 
we neglect it w.r.t. Oixoo). The JC n which are not written explicitly will give negligible contribution, see below. 
It is convenient to introduce new variables a and b such that 

ti =t{a + l/4+b), t 2 =t(a+l/A-b), (12) 

then 

Xoo = t (a - b 2 ), dtidt 2 = 2t 2 dadb, (13) 
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the condition Xoo > is equivalent to a > 0, b 2 < a and ~ t means a,b ~ 1. 

Consider the quantities A n ^j. Taking into account ifS]). the low-i behaviour of A n> ij, for all three amplitudes, is 



(14) 



Comparing with (jlip . we conclude that A\K\, A2K2 and A4_/C4_ give a contribution of order 0{y/—i) w.r.t A3/C3 
and A4+/C4+, so further we will neglect the former and consider the latter. The explicit expressions for the leading 
terms in A n ^j are, 
for the amplitude Gm- 
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for the amplitude Qm'- 
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and for Qe- 
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where we use matrix notation A ri 
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for compactness. 



Now let us consider the FFs Fi(ti)Fj(t 2 ) = Fi(ta+t/4+tb)Fj(ta+t/4-tb). For a > 1 we have b < -Ja < a, so 
Fi{ti)Fj(t2) ~ Fi(ta + t/4)Fj(ta + t/4). In the main region of interest a~ 1 wc expand F in Taylor series around 
ta + t/4 and obtain 



2F i F J 



O (t/t 



FiitjFjfa) - F^Fjih) « i6(i^ - F^) + O (tAo) 3 , 
where FFs and their derivatives in the r.h.s are taken at ta + i/4, and t is characteristic scale for FFs, 

to = 6/(r 2 ) w 0.3 GcV 2 for proton. 



(18) 
(19) 

(20) 



Note that to "C 4M 2 (this is because proton radius is related to pion rather than proton mass), and thus an additional 
requirement arises: (t/to) 2 <C 1. 

The r.h.s. of Eq. (fP9|) is of order t/to- In the presence of FF scaling, F 2 /F\ as const, it becomes much smaller, since 
F!Fj — F'-Fi = Fi Fj (In Fi / Fj)' , so we will assume (fT9"f to be negligible and thus put everywhere 



Fi{ti)F 3 {t2) = Fi(ta + t/A)F j (ta + t/4). 



(21) 
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Replacing F(ta + t/A) by just F(0) would be too rough an approximation, since the FFs suppress the integrand at 
a > 1. 

After all the approximation made the FFs do not depend on b and the integration over b can be done analytically: 
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a - 1/4 + V>* 2 - ly/a 



Collecting all the expressions together, we obtain finally 
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SGm = — x 



SGe = —x 



[J 3 + (a - l/4)(a + 3/4) J 4 - J' A ) \ F X {F X + F 2 )da, 
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(24) 

(25) 
(26) 
(27) 



where Fi >2 = Fi i2 (ta + 

The approximate formulae (|25H27p should work well for 



t/Am 2 > 1, V~t/2M < 1, (t/t Q ) 2 < 1. 



(28) 



The infrared problems would arise at t\ = 0, £2 = t or <i = t, t 2 = 0, that is, at a = 1/4. However the singularity 
in (|25H27j) at a = 1/4 is integrable, so the IR divergent terms do not appear. Actually such terms have the form 
In In A 2 w In A 2 and are negligible in our approximation. 

Consider the limit t = 0. In this case F\ = 1, F 2 = /a — 1, where is magnetic moment. In other words, the result 
will be the same as for point-like particle. Taking into account that x > 0, we perform the integration and obtain for 
t = 



~ x+l/2 x+l/2 
0(jx M = MOTT — , dyM — l-ian 5 
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(29) 



The last result for SQe corresponds to well-known second-order correction to Coulomb scattering Q. 
One point is to be clarified here. The TPE contributions should be odd functions of x From this the authors 
of Ref.[7J make a (erroneous) conclusion, that TPE amplitudes should have the form 



sg 



.2/.- 



(30) 



fc=0 



which is apparently not the case for (|29| . The source of error is that the series (|30)) converges only in certain circle 
around x = in the complex x plane. The radius of the circle is the distance to the nearest singularity of 5Q{x). 
Such singularities are branching points at s = (m + M) 2 and at u = (m + M) 2 , that is, at v = s — u = zk(AMm + i); 

neglecting m, this yields ±x ~ \J ik[ z__ t < 1- But the physical values arc x > 1, sec Eq.©. They therefore lie outside 

the convergence circle, and Eq. (|3"0"l) is not valid for these x. 

The results of numerical evaluation of TPE amplitudes according to Eqs. (|25H2"T|) arc shown in Fig. [5] (the amplitudes 
SGm and SGm are divided by fi = 2.79). The calculation is done with the dipole FFs. To check the sensitivity to the 
FF parameterization, we repeated the calculation using FF fits @; the resulting changes in SQ were small. 

The TPE amplitudes are falling rapidly near t — 0; especially this pertains to the "magnetic" corrections SGm and 
SGm- Since the determination of rms radius involves FF derivative rather than FF itself, our results suggest that 
magnetic TPE correction may be important here. Roughly speaking, TPE corrections are to be subtracted from the 
experimentally measured FFs, resulting in slower decrease of "corrected" FFs. So wc expect that inclusion of TPE 
should lessen proton rms radii. Surprising ly, Refs.[i]i, where TPE effects were calculated merely numerically, claim 
the increase of charge radius, by (0.008 - 0.013) fm [J and by 0.0015 fm 0]. 



5 




o r 1 

0.02 0.04 0.06 0.08 0.1 



Q 2 , GeV 2 

FIG. 2: TPE amplitudes at e = 0.5 (a) and at Q 2 = 0.05 
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2 (b); 5Qe (solid), SQm/h (dashed) and SGm/fj, (dotted). 
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FIG. 3: TPE corrections to the cross-section at e = 0.5, electric (dash-dotted), magnetic (dashed) and total (solid). 



In Fig. [3] we show the TPE corrections to the cross-section at e = 0.5, electric (coming from 5Qe), magnetic (coming 
from SQm) and total. We see that magnetic contribution is of the same order as the electric one and represents an 
essential piece of the total correction. At lower e it becomes even more important. 

The formalism developed here can be applied to electron- neutron scattering as well. In this case Ge,Fi ~ t, and 
we see from (|25I27[) that SQm ~ t and SGe ~ t 2 . In other words, the TPE amplitudes are suppressed, with respect 
to the Born amplitude, by a factor of at, which is much smaller than just a. Consequently, TPE corrections to the 
elastic en scattering are negligibly small at low t. 

In summary, we obtained simple approximate formulae for TPE amplitudes of elastic electron-proton scattering 
at low Q 2 . Numerical calculations indicate that the magnetic TPE corrections are relatively large and have a sharp 
Q-dcpendcncc, which can affect the determination of proton rms radii. 
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